Mon. Not. R. Astron. Soc. 000, [TV?? (2002) Printed 17 March 2010 (MN WTtfi. style file v2.2) 



Simulations on High- z Long Gamma-Ray Burst Rate 
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ABSTRACT 

Since the launch of Swift satellite, the detections of high-z (z > 4) long gamma-ray 
bursts (LGRBs) have been rapidly growing, even approaching the very early Universe 
(the record holder currently is z—8.3). The observed high-z LGRB rate shows sig- 
nificant excess over that estimated from the star formation history. We investigate 
what may be responsible for this high productivity of GRBs at high-z through Monte 
Carlo simulations, with effective Swif/HAT trigger and redshift detection probabilities 
based on current Swift/BAT sample and CGi?0/BATSE LGRB sample. We compare 
our simulations to the Swift observations via logiV — logP, peak luminosity (L) and 
redshift distributions. In the case that LGRB rate is purely proportional to the star 
formation rate (SFR), our simulations poorly reproduce the LGRB rate at z > 4, al- 
though the simulated log N — log P distribution is in good agreement with the observed 
one. Assuming that the excess of high-z GRB rate is due to the cosmic metallicity 
evolution or unknown LGRB rate increase parameterized as (1 + z) A , we find that 
although the two scenarios alone can improve the consistency between our simula- 
tions and observations, incorporation of them gives much better consistency. We get 
0.2 < e < 0.6 and 5 < 0.6, where e is the metallicity threshold for the production of 
LGRBs. The best consistency is obtained with a parameter set (e, <5)=(~ 0.4, ~ 0.4), 
and BAT might trigger a few LGRBs at z ~ 14. With increasing detections of GRBs 
at z > 4 (~ 15% of GRBs in current Swift LGRB sample based on our simulations), 
a window for very early Universe is opening by Swift and up-coming SVOM missions. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are the most luminous events 
in the Universe. They may be detected out to distance of 
z > 10 (Lamb & Reichart 2000; Bromm & Leob 2002; 
Gou et al. 2004; Lin et al. 2004). It is believed that long 
GRBs (LGRBs) with T 90 > 2 s (Kouveliotou et al. 1993) 
are from the death of massive stars (e.g., Woosley 1993; 
Paczynski 1998; Woosley & Bloom 2006), hence their rate 
may trace the star formation rate (SFR) (Totani 1997; Wi- 
jers et al. 1998; Lamb & Reichart 2000; Blain & Natara- 
jan 2000; Porciani & Madau 2001; Zhang & Meszaros 2004; 
Piran 2004; Zhang 2007). Simulations for this scenario can 
roughly reproduce the observed luminosity (L), redshift, and 
log N — log P distributions (e.g., Liang et al. 2007; Virgil et 
al. 2009a, b), however, evidence that LGRBs do not follow 
star formation unbiasedly has been growing with rapidly in- 
creasing number of high-z GRBs detected with Swift/BAT 
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(e.g., Daigne et al. 2006; Le & Dermer 2007; Yiiksel & Kistler 
2007; Salvaterra & Chincarini 2007; Guetta & Piran 2007; 
Salvaterra et al. 2008; Kistler et al. 2008, 2009; Wang & 
Dai 2009). After the detection of GRB 050904 at z = 6.29 
(Kawai et al. 2006; Haislip et al. 2006) three year later, other 
two higher-z LGRBs were detected, GRB 080913 at z ~ 6.7 
(Greiner et al. 2009) and GRB 090423 at z ~ 8.3 (Tanvir et 
al. 2009). This ever-increasing number of detectable high-z 
LGRBs shows that the LGRB rate at high redshift may be 
higher than that expected previously. Kistler et al. (2008) 
found significant excess of the LGRB rate over the SFR at 
z ~ 4 based on a sample of 36 luminous LGRBs. They sug- 
gested that some unknown mechanisms lead to this observed 
enhance in the GRB rate, such as the cosmic evolutions of 
the metallicity history and LGRB luminosity function. 

Cosmic chemical evolution predicts that LGRBs should 
occur preferentially in metal poor environment. From a the- 
oretical point of view, this is not surprising since lower 
metallicity leads to weaker stellar winds and hence less an- 
gular momentum loss, more probably resulting in rapidly 
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rotating cores in the stars, as required in the collapsar 
model for GRBs (e.g., Woosley 1993; MacFayden & Woosley 
1999; Yoon & Langer 2005). With host galaxy observations 
(Fynbo et al. 2003; Conselice et al. 2005; Gorosabel et al. 
2005; Chen et al. 2005; Starling et al. 2005; Fruchter et al. 
2006; Fynbo et al. 2006; Mirabal et al. 2006), it is found that 
the metallicity of LGRB hosts are on average lower than that 
of core-collapse supernovae (Wolf & Podsiadlowski 2007). 
The Swift sample also suggests a modest propensity for low- 
metallicity, evidenced by an increase in the rate density with 
redshift (Butler et al. 2009) . One major restriction of LGRBs 
to low metallicity is that LGRBs do not follow star forma- 
tion in an unbiased manner and LGRB rate may peak at 
a significantly higher redshift than supernovae (Lin et al. 
2004; Langer & Norman 2006; Firmani et al. 2006). Assum- 
ing that LGRB rate traces both the SFR and the metallicity 
evolution history, Li (2008) found that the observed redshift 
distribution of the Swtft LGRBs can be reproduced with a 
fairly good accuracy. 

It was also suggested that LGRBs may have experienced 
some sort of luminosity evolution with redshift, being more 
luminous in the past (Lloyd-Ronning et al. 2002; Firmani 
et al. 2004; Kocevski& Liang 2006; Salvaterra et al. 2008; 
Salvaterra et al. 2009) or having a higher efficiency of LGRB 
production rate by massive stars at high-z (Daigne, et al. 
2006) . This effect may also modify the observed LGRB rate. 

This paper is dedicated to exploring what may be re- 
sponsible for the high productivity of LGRBs at high-z 
through simulations. Four scenarios are tested: (A) A LGRB 
rate that purely follows SFR; (B) a LGRB rate that follows 
SFR in cooperation with some sort of LGRB rate evolv- 
ing with redshift characterized as (1 + z) s (Kistler et al. 
2009); (C) a LGRB rate that follows SFR in cooperation 
with cosmic metallicity evolution (taking the form as Langer 
and Norman 2006); (D) a LGRB rate that follows the SFR 
in conjunction with both LGRB rate evolution and cosmic 
metallicity evolution. Throughout this paper a flat Universe 
with H = 71 km s _1 Mpc _1 , Q M =0.3 and Q A =0-7 is as- 
sumed. 



2 MODEL 

The number of LGRBs per unit time at redshift z ~ z + dz 
with luminosity L ~ L + dL is given by 



during an observational period T should be 



dN 
dtdzdL 



Rlgrb(z) dV{z) 



l + z 



dz 



<£>(L), 



(1) 



where 7?lgrb(«) is the event rate of LGRB (in units of 
Gpc~ 3 yr _1 ) as a function of z, (1 + z) _1 accounts for the 
cosmological time dilation, $(1/) is the LGRB luminosity 
function, and 



dV_ 

dz 



Ho (l + z)2[fi M (l + z)3 + fi A ] 1 /2 



(2) 



is the co-moving volume element at redshift z in a flat A 
CDM universe. The Dl(z) is the luminosity distance at z. 
Considering an instrument with flux threshold P tb and an 
average solid angle fi for the aperture flux, the number of 
LGRB triggers with peak luminosity in a range [L m i n , L m „] 



N = 



4tv 



JL min JO 



Rlgrb(z) dV{z) 



l + z 



dz 



dz, 
(3) 

where r)t(P) is the trigger probability of a burst with ob- 
served peak flux P, given by P — L/A-KD\(z)k, where 
the k factor corrects the bolometric flux in the burst rest 
frame (1 — 10 4 keV in this analysis) into the observed flux 
in an instrument energy band with the spectral informa- 
tion. The z max of a given burst with luminosity L is de- 
termined by the instrumental flux threshold P t h through 
Pth = L/47rDi(z max )fc. 

As discussed in §1, we express the LGRB rate as a func- 
tion of redshift with 



Plgrb(z) oc RsFn(z)f(z)Q(e,z) 



(4) 



where Rsfr(z) is the star formation rate, /(z) accounts 
for potential cosmological evolution effects, 6(e, z) is the 
fractional mass density belonging to metallicity below zZq 
at a given z (Zq is the solar metal abundance), and e is 
determined by the metallicity threshold for the produc- 
tion of LGRBs. In our simulations, /(z) is simplified to 
f(z) = (1 + z) s as adopted by Kistler et al. (2009), and 
the 0(e, z) follows the law (Langer & Norman 2006), 



6(e,; 



(5) 



where k = —1.16 is the power-law index in the Schechter 
distribution function of galaxy stellar masses (Panter et al. 
2004), /3 = 2 is the slope in the linear bisector fit to the 
galaxy stellar mass- metallicity relation (Savaglio 2006), and 
T(a,x) and T(x) are the incomplete and complete gamma 
function, respectively. The value of e varies from 0.2 to 0.6 
(Modjaz et al. 2008, Li 2008). 

The star formation history is not well constrained by 
data at high redshifts (z > 4). The current observed SFR 
rapidly increases at z < 1, keeps almost a constant up to 
z = 4.0, then drops at z > 4. The sharp drop at z > 4 
may be due to large dust extinction at such high redshifts. 
We adopt a SFR parameterized as (Rowan-Robinson 1999; 
Hopkins & Beacom 2006), 

R M-B j( 1 + z ) 3A4 z^z peak , 
Rsfr(z) - Ho < 3 44 (6) 

I \i r Zp e ak ) Z S> Zp ea k> 

where z pea k = 1 and Ro is a normalization parameter. The 
slope of the segment z < z pea k is the well-fitted parameter 
in Hopkins & Beacom (2006). 

It was suggested that low luminosity GRBs (LL-GRBs) 
may be a distinct GRB population from high-luminosity 
GRBs (HL-GRBs) with much higher local event rate (Soder- 
berg et al. 2004; Cobb et al. 2006; Liang et al. 2007; Chap- 
man et al. 2007). Liang et al. (2007) proposed a two- 
component broken power-law to model the GRB luminosity 
function. They characterize the luminosity function for each 
GRB population as 



(7) 



where <f>o is a normalization constant to assure the integral 
over the luminosity function being equal to unity, and L b is 
the break luminosity. We do not consider the cosmological 
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evolution of the luminosity functiorQ. The global luminosity 
function thus is given by 

<f(L) = <f [Po L ^ LL (L) + P% L ^{L)] (8) 

where $o is a normalization constant, pg L an d Po L are the 
local rates of low and high luminosity LGRBs, respectively. 



3 SWIFT GRB SAMPLE AND SWIFT/BAT 
TRIGGER PROBABILITY 

In order to eliminate the selection effect of different in- 
struments, we use only the LGRBs detected by Swift/BAT 
to form a homogeneous sample. As of 3f May 2009, Swift 
had detected 385 LGRBs, with redshift measurement for 
124 GRBs. The peak fluxes, measurement redshift and the 
spectral information for the GRB in our sample are taken 
from the NASA websitfl The peak fluxes are measured on 
1-s timescale. It is known that the GRB spectrum is ad- 
equately fit with a smooth broken power-law, the so-called 
Band function (Band et al. 1993). Because of the narrowness 
of the BAT energy band, the BAT data cannot adequately 
constrain the spectral parameters of GRBs, and the GRB 
spectra observed with BAT are well fitted with a simple 
power law, N tx v~ T (Zhang et al. 2007). Empirically, F 
is roughly correlated to the observed peak spectral energy 
E p by logSp = (2.76 ± 0.07) - (3.61 ± 0.26) log T (Zhang 
et al. 2007; Sakamoto et al. 2009). We employ this empiri- 
cal relation to estimate the E p of Swift GRBs and correct 
the observed peak luminosity to a bolometric band (1 — 10 4 
keV) assuming the low and high energy spectral indexes as 
Ti = 0.83 and T 2 = 2.35 (Preece et al. 2000; Kaneko et al. 
2006) for all Swift GRBs and simulated bursts. 

Selection effects involved in a GRB sample may disguise 
its truly intrinsic distribution (Bloom 2003; Fiore et al. 2007; 
Kistler et al. 2008; Le & Dermer 2007; Coward 2007; Coward 
et al. 2008). Two kinds of selection effects are concerned in 
our simulations. One is the Swift/BAT GRB trigger and the 
other is the redshift determination through spectroscopy of 
the optical/NIR afterglow or of the GRB host galaxy. The 
first one affects the peak flux distribution of a GRB sample 
and the other one affects the L — z distribution. 

The Swift/BAT trigger is quite complex and its sensitiv- 
ity for GRBs is very difficult to parameterize (Band 2006). 
Practically, a GRB with peak flux slightly over the instru- 
ment threshold would not always make a trigger. The BAT 
trigger rate is about 64 LGRBs per year per solid angle, 
roughly equal to that of BATSE (~ 67 LGRBs per year 
per solid angle, see Stern et al. 2001). Comparison of the 
peak flux distributions of the LGRBs detected with BAT 
and BATSE is shown in Fig. [T] The Kolmogorov-Smirnov 
Test (K-S test) yields a probability pk-s = 0.79, indicating 
that the two distributions are statistically indistinguishable. 



1 The luminosity evolution with being luminous in the past 
(Lloyd-Ronning et al. 2002; Firmani et al. 2004; Koccvski& Liang 
2006; Salvaterra et al. 2008; Salvaterra et al. 2009) may also re- 
sult in an increase of GRB rate with redshift as (1 + z) s (Kistler 
et al. 2008). 

2 http: / /swift . gsfc . nasa.gov /docs/ swift /archive /grb.table / 



Therefore the two detectors should have comparable sen- 
sitivity for LGRBs. We consequently use the same trigger 
probability of Swift/BAT as that of CGPO/BATSE. 

Stern et al. (2001) scanned the off-line daily archival 
data from BATSE and identified 3713 LGRBs, of which 1916 
LGRBs are triggered events. Assuming that the peak fluxes 
of both trigged and non-trigged GRBs are over the BATSE 
threshold, we define the trigger efficiency rjt of a GRB as a 
function of P with the ratio of trigger event number to the 
total detected event number (including both triggered and 
non-trigged events) in a peak flux bin P + dP. The result is 
shown in Fig. [2] where P is converted to Swift energy band 
[15, 150] keV with the observed spectral parameters. It is 
found that the instrument trigger threshold is roughly 0.2 
photon s _1 cm -2 . r\ t as a function of P can be parameterized 
as 

(5.0P 3 85 ,P < 0.45 
Vt ~ [0.67(1.0 -0.40/P) 052 ,P > 0.45 ' l<) ' 

Note that r/t remains almost a constant at 67% for P > 1 
counts cm -2 s _1 . Some non-triggered events are very bright 
(the strongest one is 24 ph cm ~ 2 s _1 ). This is due to the 
dead time when the trigger was disabled during data read- 
outs or when CGRO passed through regions of very high 
ionospheric activity (Stern et al. 2001). This effect does not 
affect our our simulations. We use this trigger probability 
curve for Swift/BAT in our simulations. 

Fiore et al. (2007) pointed out that about 30% of Swift 
GRBs have measured redshifts (124/385 in the current Swift 
GRB sample). In order to simulate a GRB sample with 
redshift measurement, we empirically model probability of 
redshift measurement for a BAT triggered burst. From the 
current GRB sample, redshifts are preferentially measured 
for brighter GRBs at lower redshifts. We show the redshift 
measurement probability r\ z as a function of log(P) in Fig. 
[3] where the probabilities are calculated by the ratio of the 
number of GRBs with redshift measurement in each bin to 
the number of triggered GRBs in the corresponding bin. It 
is found that the probability of redshift measurement for 
the burst in our sample does not strongly depend on P. In 
fact, the redshift measurement is complicated, depending on 
many artificial effects, such as the optical follow-up, spectral 
line detection, etc (Bloom 2003). Even though, we still use 
an empirical function parameterized as 

r) z = 0.26 + 0.032e 1611ogP (10) 

to simulate GRB samples with redshift measurement. 



4 MONTE CARLO SIMULATIONS 

Based on the model discussed above, we make simulations 
to re-produce the Swift observations. In our simulations, 
the model parameters are constrained with comparison of 
our simulations to observations. Note that the parameters 
of the LL-GRB luminosity function are similar to that re- 
ported by Liang et al. (2007), i.e., a\ L = 0, a£ L = 3.0, 
and = 7.5 x 10 46 erg s" 1 . For HL-GRBs, we find that 
Lf L = 2.75 x 10 52 erg s~ x and af L = 1.36 can well repro- 
duce the observations at L < L^ L . This is roughly consistent 
with that reported by Stern et al. (2002), Schmidt (2001), 
and Lloyd-Ronning et al. (2004). We fix these parameters 
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in all cases we study below. The parameter a^ L are slightly 
adjustable, varying from 2.0 to 2.5, in order to get the best 
consistency between our simulations and observations for 
different cases. The Pq L is ~ 1 Gpc -3 year -1 (e.g., Schmidt 
2001; Stern et al. 2002; Lloyd-Ronning et al. 2004; Liang 
et al. 2007), but it was suggested that po L = 100 ~ 1000 
Gpc -3 year -1 , much higher than p^ L (e.g., Cobb et al. 2006; 
Liang et al. 2007; Chapman et al. 2007). The pg L is very 
uncertain, but it does not significantly affect the K-S test 
results since only a few low-luminosity GRBs in the simu- 
lated GRB sample. The e and 8 that describe the metallicity 
history and the GRB rate evolution, respectively, are free 
parameters. We search for the best model parameter sets by 
evaluating the consistency between the simulated P, L, and 
z distributions and the observed ones with the K-S test. A 
larger value of pk-s indicates a better consistency. A value 
of Pk-s > 0.1 is generally acceptable to claim the statisti- 
cal consistency, and a value of Pk-s < 10 -4 convincingly 
rejects the hypothesis of the consistency (e.g. Bloom 2003). 
For simulations for a given set of parameters, we first evalu- 
ate the consistency of the log TV — log P distribution for the 
simulated triggered sample since the L and z distributions 
suffer from the biases of redshift detection. By screening the 
simulated triggered sample with the redshift measurement 
probability curve the consistencies of the P, L and z distri- 
butions should be significantly improved (with pk-s > 0.1) 
for a reasonable parameter set. The details of our procedure 
are described as follows. 

(i) Simulate a GRB with luminosity L at redshift z, 
GRB(L,z). The redshift is generated from the co-moving 
number density of GRBs at redshift z + dz, i.e., 

-Rlgrb(z) dV 



H(z) = 



1 + z dz 



(11) 



and L is simulated with the probability distributions derived 
from Eq.[i We take [L min , L max ] = [10 45 , 10 55 ] erg s" 1 and 
a redshift range of z = — 20. 

(ii) Derive E p in the burst frame with the L — E p relation, 



E p » 200keVL^ 2 /C. 



(12) 



where C is taken from a random distribution in the range 
of [0.1, l](Liang et al. 2004). The photon indices prior and 
post the break energy are assumed to be 0.83 and 2.35, re- 
spectively. With the spectral information, we calculate the 
peak photon flux P in the BAT band that corresponds to an 
energy band [15 x (1 + z), 150 x (1 + z)] in the burst frame. 
The simulated GRB then is characterized as GRB (L, z, P). 

(iii) Screen a mock GRB with the BAT trigger probability 
(Eq. [9]l . We get the trigger probability r\ t of burst with peak 
flux P from Eq. [9] and generate a random number Qt in 
the range (0,1). If < Qt < r/ t , we pick up this event as a 
triggered GRB. 

(iv) Repeat the above steps to simulate a BAT-triggered 
GRB sample of 385 GRBs and evaluate the consistencies of 
its log TV — log P distribution with the BAT observations by 
the K-S test. 

(v) Screen the simulated triggered sample (385 GRBs) 
with the redshift measurement probability curve (Eq. I10[) . 
To do so, we calculate the r\ z value for a given burst, and 
generate a random number Q z in the range (0,1). If < 
Qz < r] z , we pick up this event as a triggered GRB with 
redshift detection. 



(vi) Evaluate the consistencies of log TV — logP, L and 
z distributions of the simulated known-redshift by the K-S 
test. In a reasonable parameter set these distributions should 
be greater agreement with the observations than that of the 
simulated triggered sample, i.e., with larger pk-s values. 

We create simulations for four cases as below. We 
present our selected simulation results in the following sec- 
tion. The selected model parameters and K-S test results are 
summarized in Table 1, and distributions are shown in Figs. 
4-6 and 8. Our simulations show that the highest redshift of 
LGRBs that may be triggered with Swift/B AT is ~ 14. 



4.1 Case A: Rlgrb(z) oc Rsfh(z) 

This case is for the GRB rate that purely follows the SFR. 
This was performed by Virgili et al. (2009a). We make im- 
provements in this case by considering the trigger proba- 
bility and redshift detection probability in our simulations. 
Comparisons between simulations and Swift observations are 
shown in Fig. [4] It is found that our simulations well re- 
produce the observed log TV — log P distribution, yielding 
Pk-s = 0.423 and pk-s = 0.585 for the trigged and known- 
redshift samples, respectively. However, the observed lumi- 
nosity and redshift distributions are poorly reproduced. The 
simulated GRBs at z ~ 1 is significantly over-produced, 
but the GRBs at z > 2 have significant deficit. This case 
also over-produces the GRBs around L = 10 erg s - , and 
under-generates GRBs with L > 10 53 erg s - . The px-s for 
the logL and logz distributions are smaller than 0.01. 



4.2 Case B: Rlgrb(z) oc 7?s FR («)e(e, z) 

This case assumes that the GRB rate is proportional to the 
star formation history incorporating with the cosmic metal- 
licity history. We find that e ~ 0.5 gives the best consistency. 
Our results are shown in Fig. [5] In general, the simulations 
for this case are more consistent with the observations than 
Case A. The log TV — log P distribution is well reproduced, 
with pk-s = 0.372 and pk-s = 0.449 for the simulations on 
the trigged and known-redshift samples, respectively. The 
consistency of the simulated luminosity distribution is also 
accepted, withp_R-_s = 0.052 andp^-s = 0.292 for the trig- 
gered and the known-redshift samples, respectively. How- 
ever, the simulated z distributions are only marginally ac- 
ceptable, with pk-s = 0.033 and pk-s = 0.009 for the sim- 
ulations on the trigged and known-redshift samples, respec- 
tively. Note that the consistency of z distribution becomes 
much worse for the simulated known-redshift sample. 



4.3 Case C: Rlgbb(z) oc R S Fn(z)f(z) 

This case assumes that the GRB rate follows the SFR in 
conjunction with some sort of GRB rate evolving with red- 
shift characterized by (1 + z) s (Kistler et al. 2009). We 
find that S — 0.85 can yield the best consistency. The re- 
sults are shown in Fig. [6] This case well reproduces the ob- 
served log TV — logP, z, and L distributions, with pk-s = 
0.889, 0.242, 0.369, respectively, for the simulations on the 
known-redshift samples, respectively. However, the consis- 
tency of the z and L distributions of the simulations for 
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the triggered sample are only marginally acceptable (with 
Pk-s < 0.1). 

4.4 Case D: Rlgrb(z) oc Rsfr(z) f(z)<d(e, z) 

In this case we assume that the GRB rate follows the SFR 
incorporating with both GRB rate evolution and cosmic 
metallicity evolution. We search for the best 8 and e by 
means of a global K-S test probability defined as Pk-s = 
Pk-s x Pk-s x Pk-s, where p£_ s , p K -s, and p h K _ s are 
the probabilities of the K-S tests for the logiV — logP, z, 
and L distributions of the simulated GRB sample with red- 
shift measurement. The contours of pk-s J? 0.1 are shown 
in Fig. □ We find that 0.2 < e < 0.6 and 8 < 0.6 in case 
of pk-s > 0.1. The parameter set (e, 8)=(~ 0.4, ~ 0.4) 
produces the best consistency between our simulations and 
observations, as show in Fig. \E\ yielding Pk-s = 0.698, 
Pk-s — 0.294 and Pk-s = 0.209 for the triggered sample 
and Pk-s = 0.993, Pk-s = 0.843 and p K _ s = 0.966 for the 
known-redshift sample. 



5 DISCUSSION 

As shown above, the observed excesses of both high-z GRB 
rate at z > 4 and luminous GRBs with L > 10 53 erg s _1 over 
that predicted by assuming that the GRB rate simply traces 
the star formation rate (Case A), can be explained with the 
cosmic metallicity evolution (Case B) and an unknown GRB 
rate increase as (1 + z) s (Case C). 

Our simulations of Case D that incorporates Cases B 
and C well reproduce the observations. Taking this case as 
a preferred one, we find that the percentage of GRBs with 
2 > 4 in the current BAT GRB sample is ~ 15% (~ 58 
bursts) and ~ 4% of the GRBs (~ 15 bursts) have z > 6, 
consistent with that reported by Butler et al. (2009), but 
slightly larger than that predicted by (Daigne, et al. 2006). 
Swift has operated for five years. We thus simulate the obser- 
vations for a 10 year operation. A uniform sample of ~ 1000 
GRBs with ~ 300 redshift detections (in the current follow- 
up status) would be established in this mission period. Our 
simulations show that Siui/t/BAT-like instruments can trig- 
ger GRBs at redshift up to 14, and ~ 0.5% triggered Swift 
GRB may be at redshift z > 10. Exciting evidence is from 
Swift GRB 100205. The detected high red color of this GRB 
in its afterglows, K-H=1.6±0.5 mag, likely suggests that its 
redshift is 11 z sC 13.5, if this red color is due to Lyman- 
q absorption within the H filter bandpass (Cucchiara et al. 
2010). The era of exploring the deepest Universe with GRBs 
is coming, e.g., with rapid follow-up and localization capa- 
bilities on the ground and the space-based multi-band as- 
tronomical variable object monitor (SVOM) mission that 
is being developed in cooperation between the China and 
French (Gotz et al. 2009; Basa, et al. 2008). Thanks to the 
low energy trigger threshold (~ 4 keV) of the ECLAIRs on- 
board, SVOM is more sensitive than previous missions for 
the detection of soft, hence potentially most distant, GRBs. 
The on-board visible telescope (VT) would quickly refine 
the GRB position for three ground based dedicated instru- 
ments, two robotic telescopes (GFTs) and a wide angle op- 
tical monitor (GWAC). Such optimized observing strategy 
would increase the redshift detection rate for high-z GRBs. 



6 CONCLUSIONS 

Motivated by the increasing detections of high-z GRBs, we 
have investigated what may be responsible for the high pro- 
ductivity of GRBs at high-z through Monte Carlo simula- 
tions. We parameterize the trigger probability of Swift/BAT 
using a large GRB sample observed with CRGO/BATSE. 
The redshift measurement probability is also derived from 
current Swift GRB sample. We compare our simulation re- 
sults to the Swift observations with log N — log P and L — z 
distributions. We show in the case that the GRB rate is 
purely proportional to the SFR, our simulations poorly re- 
produce the GRBs rate at z > 4. We explain this GRB 
rate excess over that predicted by the SFR with the cosmic 
metallicity evolution effect and unknown GRB rate increase 
as (1 + z) s . Although the two scenarios can make better 
consistency between our simulations and observations, ei- 
ther one cannot simultaneously reproduce the observations 
alone. Incorporation of the two scenarios gives greater agree- 
ment between our simulations and observations, indicating 
that the combination of both GRB rate evolution and cos- 
mic metallicity evolution would result in the observed high-z 
GRB rate excess over that predicted from the star forma- 
tion history. We get 0.2 < e < 0.6 and 8 < 0.6 in case 
of Pk-S > 0.1. The parameter set (e, 8)=(~ 0.4, ~ 0.4) 
produces the best consistency between our simulations and 
observations. Our simulations show that the percentage of 
GRBs with z > 4 in the current BAT GRB sample is ~ 15%, 
i.e., ~ 58 bursts, and a few triggered LGRBs may be at 
z ~ 14. This is sufficient to make GRBs as promising probe 
for the high- z Universe. 
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Figure 1. Comparison of the peak flux distribution of LGRBs 
detected by BATSE (dashed line) to BAT GRB sample (solid 
line). BATSE data are from Stern et al. (2001) and BAT data are 
from NASA website. The peak fluxes of the BATSE GRBs are 
converted to the BAT energy band of [15, 150] keV for compari- 
son, with a typical spectral parameter set of Ti = 0.83, T2 = 2.35 
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Figure 3. Redshift measured probability as a function of the 
peak flux with our best fit (dotted curve) for BAT LGRBs. The 
best fit curve is parameterized in Eg. I10I 



1 10 100 

P (photons cm" 2 s" 1 ) 



Figure 2. Trigger probability as a function of the peak flux with 
our best fit (dotted curve) for BATSE LGRBs. The best fit curve 
is parameterized in Eq. [9] 
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Table 1. The selected model parameters from our simulations and the corresponding K-S test probability. 



Case Model parameters" K-S test 6 





e 


<5 






Pk-S 


Pk-S 


Pk-S 


A 






50 


2.5 


0.423 | 0.585 


< 0.001 | < 0.001 


< 0.001 | 0.005 


B 


0.5 




50 


2.2 


0.372 | 0.449 


0.033 | 0.009 


0.052 | 0.292 


C 




0.85 


300 


2.0 


0.258 0.889 


0.060 | 0.242 


0.018 0.369 


D 


0.1 


0.4 


200 


2.0 


0.698 0.993 


0.294 0.843 


0.209 | 0.966 



Note: a The parameters of the luminosity function of LL-GRBs are fixed at a\ L = and = 3.0, L^ L = 7.5 X 10 46 erg s _1 . 
For HL-GRBs, Lf L also fixed at 2.75 X 10 52 erg s" 1 and af L fixed at 1.36. b From left to right: The K-S test probabilities of 
log N — logP, log(l + z). and logL distributions (triggered GRBs | z-known sample). 
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Figure 4. Comparisons of two-dimensional log L — log z distri- 
butions and one-dimensional logL, log 2, and logP (accumula- 
tive log N — log P) distributions between the observed Swift/BAT 
GRB sample (solid dots and solid lines) and our simulations (open 
dots and dashed lines) for Case A: Plgrb( 2; ) °^ Psfr( z )- ' e /* 
four panels are for the trigged GRB sample. The contours in the 
log L — log 2 plane show the relative probability distributions of 
the simulated GRBs with different color lines: red for 99.7%, ma- 
genta for 95.5%, and orange for 68.7%. Right four panels are for 
the sample with redshift measurement. One dimensional K-S test 
probabilities for the comparisons are marked in each panel. 
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Figure 5. The same as Fig. 4, but for Case B: Rlgrb(^) oc 
R S FR(z)e{e,z). 
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Figure 6. The same as Fig. 4, but for Case C: Rlgrb (z) oc 
R S fr(z)(1 + z) s . 
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Figure 8. The same as Fig. 4, but for Case D: Rlgrb(z) 
RsFR(z)(l + z) s e(e,z). 

with a parameter set (<5, e) = (0.4, 0.4). 



Figure 7. Contours of Pk—S m the e — 8 plane for Case D: 
-Rlgrb(z) oc _R S fr(z)(1 + z) s &(e,z), where Pk-S = Pk-S x 
Pk-S x Pk-S' Px_g, pfr_«j and pj^_g are the probabilities 
of the K-S tests for the log N — log P, z , and L distributions of the 
simulated GRB sample with redshift measurement, respectively. 
The cross marks the best consistency between our simulations 
and observations, i.e., (<5, e) = (0.4, 0.4). 



